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ELECTRONS, MUONS AND HADRONS IN EXTENSIVE AIR SHOWERS 
AND HOW DO THEY DEPEND ON NUCLEAR INTERACTION MODEL 

(Part II) 

J.A.Wrotniak and G.B.Yodh 

Department of Physics and Astronomy 
University of Maryland 
College Park, MD 20742, U.S.A. 

Here we present some of the results of Monte 
Carlo simulations of extensive air showers for 
nuclear interaction models as outlined in our 
contribution HE 4.1-7 to this Conference. 

In the notation used below, numbers in brackets () denote 
mean square errors in last decimal digit units, k, M, G, T 
and E stand for appropriate powers of 10. For the scarcity 
of place, the radial data on showers are not included. 


Table 1. Average shower size at 1000 q/cm a 
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Table 2. Fluctuations of the shower size at 1000 g/cm* (s.d.of Log{base 10} Ne) 
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Table 3. Average depth of shower maximum (g/cm* ) 


1 

1 
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1 
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1 
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i 

i 
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1 
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1 
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i 
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557(7) | 
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1 
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i 
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i 
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2 P | 
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1 
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10 P | 
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1 
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639(7) 

i 
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50 P | 
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1 
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i 
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1 

729(6) 1 

1 


i 
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1 E |> 
1 
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1 
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1 

1 
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1 
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1 
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Table 4. Average shower size at maximum (exactly: geometric mean values) 


i r 

i i 
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1 

i 
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1 

"i 
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1 
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1 
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1 
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1 
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1 

1 
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1 
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I 

1 
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1 

33.3(3)M| 

1 

1 
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1 
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1 

1 
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1 1 8 ( 1 ) M 1 

1 1 
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1 
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Fig.l. The depth of maximum for some of Fig. 2. Energy per electron at the 

our models. shower maximum. 
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Table 5, Fluctuations in shower size at maximum (s.d.of Log {base 10} NeCmax]) 


i r 

i i 

1 Etev] I 


Primary protons 

1 

1 

1 

Primary iron 

1 1 
1 1 
1 1 

F-Y00 

M-Y00 1 
| 

1 

M-F00 1 

M-F01 1 
1 

R-F01 1 
1 

FF-Y00 

RM-F00 

1 1 
1 20 T I 

.10 

I 

.09 | 

1 

1 

.032 

.033 

1100 T I 

.08 

.08 | 

.08 | 

.07 i 

.07 | 

.027 

.029 

1500 T I 

.07 

i 

.06 I 

1 

1 

.021 

.022 

1 2 P 1 

.06 


.05 I 

1 

1 

.017 

.014 

1 10 P 1 

.06 

.06 ! 

.04 | 

.05 1 

.04 1 

.015 

.015 

1 50 P 1 

, .06 

1 

.04 | 

1 

1 

.014 

.011 

1200 P 1 

.07 

1 

.04 | 

1 

1 

.009 

.014 

1 1 E I 

1 1 

.07 

.04 I 

1 

.04 | 

.04 1 

1 

.02 I 

1 

.013 

.008 


Model M-F00 was also run at primary proton energy of 100 EeV (or 10 50 eV); 151 
showers at 1000 g/cm 1 and 64 at 1400 g/cm l were simulated. Number of 
electrons at 1000 g/cm 1 was 59.5(5) G, its fluctuations .04, at maximum 
(average depth 906(5) g/cm 1 ) : 66.1(5) G and .024, respectively. 


Table 6. The average muon number at 1000 q/cm 2 . 
(a) at E > 2 GeV 


E[eV] 


I 20 
1100 
I 500 
I 2 
I 10 
I 50 
1200 
I 1 


Primary protons 


i I I 1 

F-Y00 | M-YOO | M-F00 I M-F01 I R-F01 


I 


I 


,36(l)k| I . 34(l)k| I 

1 . 55 ( 3 )k I 1.53(2)k| 1.43(3)k| 1.09(2)kl 
6.6(l)k| | 6.1(l)kl I 

22.4(5)k I | 21.0(3)k| I 

90 ( 2 ) k I 9 1 ( 2 ) k I 82 ( 1 )k 1 52.9(6)kl 
3 3 7 ( 6 ) k I | 334 ( 4 ) k I | 

1.03(3)M| | 1 .08 ( 2 )M I I 

3.71(14 )M| 4 . 54 ( 9 )M| 4.14(7)M| 2.33(4)M| 

I I I I 


1 . 06 ( 2 ) k 


51.9(8)k 


2 . 38 ( 3 )M 


Primary iron j 

I 

I I 

FF-Y00 | RM-F00 I 


V 


'I 


. 64 ( 1 ) k I . 64 ( 1 )k I 

2 . 32 ( 2 ) k I 2 . 30 (2 )k I 


9 . 57 ( 6 ) k I 
34 . 2 ( 2 ) k I 


9.24(4 )k I 
32.3(3)k| 


147.6(6)kll36.1(6)kl 
61 2 ( 2 ) k I 567 ( 3 ) k I 
2 .03 ( 2 )M| 1 . 90 ( 1)M| 
7 . 87 ( 5 )M I 7.60(3)M| 


(b) at E > 200 GeV 


I' 


E[eV] 
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i 

i 


1 
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i 

i 
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i 

j. 

M- 
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■ 

i 

j. 
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J. 
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i 
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T 

1 

1 
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T 

1 


i 

2.9 


i 


i 


i 
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1 

<.l 
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T 

1 
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1 
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i 
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i 

8.5 

i 

8.6 

i 
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1 

7.5 
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T 

1 

29 


1 


i 

32 


i 


i 


i 

73 


1 

77 
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P 

1 

83 


1 


i 

97 


i 


i 


i 
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1 

223 


1 10 

P 

1 

299 


1 

330 

i 

342 


i 

254 

i 

279 

i 

692 


f 

773 


1 50 

P 

1 

1.03 

k 

1 


i 

1.24 

k 

i 


i 


i 

2.43 

k 

1 

2.71 

k 

1 200 

P 

1 

2.88 

k 

1 


i 

3.96 

k 

i 


i 


i 

7.11 

k 

1 

8.25 

k 

1 1 
1 

E 

1 

1 

10.5 

k 

1 

1 

13.6 k 

i 

i 

13.8 

k 

i 

i 

9.3 k 

i 

i 

10.7 k 

i 

i 

24.5 

k 

1 

1 

30.1 
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Primary iron 
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Table 7. The average hadron number above 2 GeV at 1000 q/cm a . 



1 

1 

1 Etev] 


Primary protons 

r 

i 

i 

Primary 

iron 

1 

1 

1 


1 

F-Y00 1 
1 

1 

M-Y00 1 
1 

1 

M-FOO I 
1 

r 

M-F01 1 
1 

R-F01 1 
1 

r 

FF-Y00 | 

1 

RM-F0Q 

1 

1 20 

T 

i 

4.3(3) | 

1 

1 

4.8(3) 1 

1 

1 

1 

2.0(2) I 

4.2(3) 

1100 

T 

34(3) I 

38(3) 1 

27(2) 1 

17(2) | 

19(2) I 

11(1) 1 

17(1) 

1 500 

T 

244(11) I 

1 

168(8) 1 

1 

1 

88(3) 1 

96(2) 

1 2 

P 

1 . 33 ( 6 ) k I 


701(32) I 

1 

1 

554(15) 1 

500(13) 

1 10 

P 

6 . 7 ( 2 ) k I 

5 . 9 ( 2 ) k I 

3 . 7 ( 2 ) k 1 

1 . 95 ( 9) k I 

1 . 65 ( 9 ) k 1 

4 . 41 ( 9 ) k I 

3 . 19 ( 6 ) k 

1 50 

P 

34 ( 2 ) k I 

1 

18 . 3 ( 7 ) k 1 

1 

1 

28 . 8 ( 4 ) k I 

18 . 0 ( 5 ) k 

1200 

P 

131 ( 4 ) k 1 

1 

65 ( 3 ) k 1 

1 

1 

131(3)k| 

78 ( 2 )k 

1 1 
1 

E 

566 ( 19 ) k 1 
1 

550 ( 15 ) k I 
1 

293 ( 12 ) k 1 
1 

130 ( 6 ) k I 
1 

106 ( 5 ) k I 
1 

711 ( 10 )k 1 
1 

388 ( 8 )k 




Fig. 4. Relative fluctuations of the shower 
size at maximum. 


[One more model was used for EAS generation: M-F10, with all 
features exactly like M-FOO, but with quite different 
multiplicity distributions. All shower characteristics 
checked by us were in statistically good agreement between 
these models ] . 

In principle, the data shown here should speak for 
themselves. We would like, however, conclude with three 
remarks: 


* The most significant part of scaling violation effect is 
generated by the inclusion of rising cross-section. 

* Among the models considered the lowest value for 
Eo/N[max] is obtained when rapidly rising cross-section 
and charge exchange are both included (model R-F01). The 
value is still 1.38 GeV/electron. 

* Except at the highest energies, the sensitivity to atomic 
mass of the primary is greater than to specific 
assumptions about multiple production. 
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